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Abstract. We investigate the microscopic model for the frustrated layered antiferromagnets 
CS2CUCI4 and Cs2CuBr4 by performing ab initio density functional theory (DFT) calculations 
and with the help of the tight-binding method. The combination of both methods provide 
the relevant interaction paths in these materials, and we estimate the corresponding exchange 
constants. We find for CS2CUCI4 that the calculated ratio of the strongest in-plane exchange 
constants J' / J between the spin-1/2 Cu ions agrees well with neutron scattering experiments. 
The microscopic model based on the derived exchange constants is tested by comparing the 
magnetic susceptibilities obtained from exact diagonalization with experimental data. The 
electronic structure differences between CS2CUCI4 and Cs2CuBr4 are also analyzed. 



1. Introduction 

The frustrated layered antiferromagnets CS2CUCI4 and Cs2CuBr4 have recently attracted a lot of 
interest because of unconventional magnetic properties [HEIE]. While in CS2CUCI4, a magnetic 
field induced Bose-Einstein condensation (EEC) of magnons [3] and spin-fractionalization [2] 
have been observed, Cs2CuBr4, though being iso-structural to CS2CUCI4, shows no sign of BEC. 
Instead, its field-dependent magnetization shows two plateaux [Ij. 

Both compounds have been classified as two-dimensional systems, with an underlying 
anisotropic triangular lattice of spin-1/2 Cu^+ ions (right panel of figured]). The ratio of the 
two antiferromagnetic exchange constants in the triangular lattice, J^/ J (see figured]) has been 
found from neutron scattering experiments to be 0.34 for CS2CUCI4 [4J and 0.74 for Cs2CuBr4 
[1] . The interlayer couplings as well as the Dzyaloshinski-Moriya interaction also present in these 
systems are suggested to be an order of magnitude smaller than J and J^ 

Yet, since extracting the exchange constants from the neutron scattering experiment data 
is based on the assumption of a suitable Hamiltonian, there is a need to support experimental 
findings with theoretical predictions. This motivated us to perform density functional theory 
calculations on CS2CUCI4 and Cs2CuBr4 as presented in this work. 



2. Computational details 

We have performed density- functional calculations in the generalized gradient approximation 
(GGA) [5j, using the full-potential linearized augmented plane wave (LAPW) method 
implemented in the WIEN2k code [6J. With RKm — 7.0 and 120 fc-points in the irreducible 
Brillouin zone, we obtained well-converged results for both CS2CUCI4 and Cs2CuBr4. For 



Cs2CuCl4 we took the experimental room-temperature structure reported in [7J, and for 
Cs2CuBr4 we used the structure reported in [8j. 

3. Crystal structure 

CS2CUCI4 and Cs2CuBr4 crystahize in the orthorhombic space group Pnma (see figure [1]). The 
unit ceh contains four formula units and has lattice parameters a — 9.769 A, 6 = 7.607 A and 
c = 12.381 A for CS2CUCI4 and a = 10.195 A, 6 = 7.965 A and c = 12.936 A for Cs2CuBr4. 
Each of the four equivalent Cu^+ ions is surrounded by four halogen (Cl~ or Br~) ions located 
at the vertices of a squeezed tetrahedron. The CuCl4~/CuBr4~ complexes arranged in a two- 
dimensional triangular lattice form layers parallel to the he plane separated along a by Cs^ ions. 
The most important superexchange pathways between Cu ions have been determined (41 Qj to 
be those within a CuCl4~/CuBr4~ layer (J and J' in the right panel of figured]). The magnetic 
interaction is mediated via halogen ions. 




Figure 1. Left: Unit 
cell of Cs2CuBr4 and 
CS2CUCI4. Right: Lattice 
of magnetic Cu sites where 
the five most important 
interaction pathways are 
shown. 



4. Electronic structure 

In panels (a) and (b) of figure [2] we present the calculated bandstructures of CS2CUCI4 and 
Cs2CuBr4. Cu t2g states that are hybridized via halogen p states lie in the range of energies 
between -0.6 eV and 0.1 eV. The twelve bands originate from the three Cu t2g orbitals which are 
split due the distorted tetrahedron crystal field of halogen ions. In CS2CUCI4, the four partially 
unoccupied higher-energy Cu ?>dxy bands are separated from the nearly degenerate Cu 'idxz and 
Cu 'idyz bands by a narrow gap. In Cs2CuBr4, the gap is absent. Out of these bandstructure 
features we conclude that the low-energy behavior of CS2CUCI4 can be described by a one-band 
model {dxy) while the three Cu t2g band model is necessary for the description of Cs2CuBr4. 




Figure 2. Bandstructure of 
CS2CUCI4 (a) and Cs2CuBr4 
(b). Dashed lines in (a) are the 
tight-binding fit. 



Both CS2CUCI4 and Cs2CuBr4 are Mott insulators. The fact that the GGA functional gives 
a metallic state is due to the insufficient treatment in this approach of electronic correlations. 



Consideration of GGA+U provides the correct insulating behavior. Since we are interested in 
obtaining the effective hopping matrix elements between the Cu ions, we will consider in what 
follows the GGA exchange-correlation functional. 

The Cu 3dxy bands of CS2CUCI4 are comparatively less dispersive than those of Cs2CuBr4. It 
is then to be expected that the hopping parameters between magnetic Cu ions in CS2CUCI4 will 
be smaller than in Cs2CuBr4. In CS2CUCI4, the bands are almost flat in the fca^-direction, which 
indicates that layers of Cu ions (parallel to the 6c-plane) are weakly coupled. Both observations 
are in agreement with experiment. 

5. Calculation of exchange couplings for CS2CUCI4 

Based on the LAPW bandstructure (figure [2] (a)), we have calculated hopping integrals between 
magnetic Cu ions of CS2CUCI4 using the tight-binding (TB) method. We constructed a single- 
orbital TB Hamiltonian intended to describe the bands near the Fermi level, which are the 
Cu 3dxy bands in CS2CUCI4. The hopping integrals t entering the Hamiltonian were varied until 
the best possible fit of the TB bands to the LAPW bands was achieved (dashed hues in figure [2] 
(a)). Application of the TB method for Cs2CuBr4 is very involved since one needs to consider 
a three-orbital Hamiltonian to describe the overlap of the three Cu t2g bands. The number of 
variational parameters of the three-band model Hamiltonian greatly increases due to all possible 
inter-orbital hoppings. 

In table [T] we list the eight hopping integrals t that we included in the Hamiltonian for 
CS2CUCI4. The interaction pathways corresponding to t, t\ ti, ^2 and ts are shown in figure 
[1] (labeled by the corresponding exchange constants J). J = + jAFM ^j^g^g jFM jg 



Table 1. Calculated hopping integrals ti with corresponding exchange constant ratios 
(Ji/Jmax) = (U/tmax.)'^ vcrsus experimentally determined exchange constants J^^^ for CS2CUCI4. 



Interaction paths (p) 


PJ 


PJ' 


PJi 


PJ2 


PJ3 


PJ4 


PJ5 


PJe 


ti (meV) 


-11.7 


-6.75 


4.20 


-5.18 


-8.39 


-2.17 


-4.11 


-3.35 


{Ji/ Juidix) ~ (^i/^max) 


1 


0.336 


0.131 


0.197 


0.516 


0.033 


0.123 


0.074 


( J / J Y^p 


1 


0.342 




0.045 











the ferromagnetic component of the exchange while J^^^ describes the antiferromagnetic 
contributions to the total exchange constant, jf^^^ can be estimated by applying perturbation 
theory on the Hubbard model as jf^^^ = ^tf/U where ti are the hopping integrals and U is the 
on-site Coulomb repulsion. An alternative way of calculating J is by considering total energies 
of various magnetic configurations. In Ref. [9J we found by considering total energy differences 
that the ferromagnetic contributions are negligible for the strongest interactions, which allows 
us to compare the obtained jf^^^ from ti with the Ji obtained from experimental data. In 
table [1] we present the calculated (second row) and experimentally fitted (third row) ratios 
Ji/Jma-x, which are independent of U. We observe that the calculated in-plane exchange ratio 
J^/J = 0.336 agrees well with the experimentally fitted 0.342. At the same time, we observe 
that the antiferromagnetic contribution of the interlayer exchange constants Ji, J2 and J3 is 
also considerable, which implies a non-negligible three dimensional contribution of the exchange 
interactions in CS2CUCI4. 

The magnetic susceptibility of a Heisenberg Hamiltonian for CS2CUCI4 based on exchange 
couplings from our DFT calculations was calculated with the exact diagonalization method (the 



Hamiltonian of a 16-site system was diagonalized with periodic boundary conditions). Two 
models were investigated: one including J and only and the other including J, J' and the 
largest interlayer coupling J3. In figure [3] we present the results together with experimental data. 
The susceptibility plot is rescaled in a way that only relative values of J's play a role [lOj. 
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Figure 3. Susceptibility of CS2CUCI4 vs. 
temperature in units of the peak temperature Tp. 
Green sohd hne: DFT-based model with J and 
J'; red dashed hne: DFT-based model with J, 
J' and J3; crosses: experimental data [TT]. The 
susceptibility is scaled as Tpx/Ag^ with A = 
NAf^%/4:k = 0.0938 cgs units and a g factor of 
^ = 2.17. 



We observe that inclusion of the interlayer coupling J3 has very little effect on the 
susceptibility and therefore other measurements like neutron scattering along the direction 
perpendicular to the be plane are needed to investigate such possible interlayer couplings. 
We note that the less than perfect agreement of our exact diagonalization calculations with 
experiment is due to finite size effects [lOj. 

In summary, we presented electronic structure calculations of CS2CUCI4 and Cs2CuBr4. 
While CS2CUCI4 can be well described by a single-band model with important J, J' intralayer 
interactions and non-negligible interlayer coupling, a three-band model is needed for Cs2CuBr4 
due to the strong involvement of all Cu t2g bands near the Fermi surface. We found that the 
magnetic susceptibility for CS2CUCI4 is not drastically changed by the inclusion of interlayer 
couplings in the model and therefore alternative measurements have to be considered to 
investigate the interlayer couplings in this material. 
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